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ABSTRACT 

I study the joint effect of dynamical friction, tidal torques and cosmological constant on clusters of galaxies formation. I show 
that within high-density environments, such as rich clusters of galaxies, both dynamical friction and tidal torques slows down 
the collapse of low-i^ peaks producing an observable variation in the time of collapse of the perturbation and, as a consequence, 
a reduction in the mass bound to the collapsed perturbation. Moreover, the delay of the collapse produces a tendency for less 
dense regions to accrete less mass, with respect to a classical spherical model, inducing a biasing of over-dense regions toward 
higher mass. I show how the threshold of collapse is modified if dynamical friction, tidal torques and a non-zero cosmological 
constant are taken into account and I use the Extended Press Schecter (EPS) approach to calculate the effects on the mass 
function. Then, I compare the numerical mass function given in Reed et al. (2003) with the theoretical mass function obtained 
in the present paper. I show that the barrier obtained in the present paper gives rise to a better description of the mass function 
evolution with respect to other previous models (Sheth & Tormen 1999 (hereafter ST), Sheth & Tormen 2002 (hereafter STl)). 
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1. Introduction 

Structure formation in universe is generated through the 
growth and collapse of primeval density perturbations 
originated from quantum fluctuations (Guth & Pi 1982; 
Hawking 1982; Starobinsky 1982; Bardeen et al. 1986 
(hereafter BBKS)) in an inflationary phase of early uni- 
verse. Density perturbations evolve towards non-linear 
regime because of gravitational instability, they break 
away from general expansion and finally they recoUapse 
at a time Tcq. The collapse time, Tco, depends on the 
characteristic of the initial fluctuation field and on the 
environment in which the perturbation is embedded. This 
last features depends on the cosmological scenario. In this 
context, the last years have shown that the most suc- 
cessful model is the ACDM model with n,„ ~ 0.3 and 
A = 0.7. The statistics of density fluctuations originated 
in the inflationary era are Gaussian and it can be ex- 
pressed entirely in terms of the power spectrum of den- 
sity fluctuations. On average the characteristics of the 
density field peaks, e.g., their mass distribution, peculiar 
velocities, etc., are completely determined by the spec- 
trum through its moments (at least during the linear and 
early non-linear phases of the collapse (BBKS)). Moreover 
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the isotropy condition imposes that all physical quanti- 
ties around density peaks is, on average, spherically sym- 
metric. However, several years ago it was realized that 
the density field distributions around the density peaks, 
which eventually will give birth to galaxies and clusters, 
depart from spherical symmetry and from the average den- 
sity profile, producing important consequences on collapse 
dynamics and formation of protostructures (Hoffman & 
Shaham 1985; Ryden 1988a,b; Heavens & Peacock 1988; 
Kashlinsky 1986, 1987; Peebles 1990). A fundamental role 
in this context is played by the joint action of tidal 
torques (coupling shells of matter which are accreted 
around a density peak and neighboring protostructures 
(Ryden 1988)), and by dynamical friction (White 1976; 
Kashlinsky, 1986, 1987, Antonuccio & Colafrancesco 1994 
(hereafter AC)). 

According to the previrialization conjecture (Peebles 
& Groth 1976, Davis & Peebles 1977, Peebles 1990), initial 
asphericities and tidal interactions between neighboring 
density fluctuations induce significant non-radial motions 
which oppose the collapse. This means that virialized 
clumps form later, with respect to the predictions of the 
linear perturbation theory or the spherical collapse model 
(hereafter SM), and that the initial density contrast, 
needed to obtain a given final density contrast, must 
be larger than that for an isolated spherical fluctuation. 
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This kind of conclusion was supported by Barrow & Silk 
(1981), Szalay & Silk (1983), Villumscn k Davis (1986), 
Bond & Myers (1993a,b) and Lokas et al. (1996). 
In particular Barrow & Silk (1981) and Szalay & Silk 
(1983) pointed out that non-radial motions would slow 
the rate of growth of the density contrast by lowering the 
peculiar velocity and suppress collapse once the system 
detaches from general expansion. Villumscn & Davis 
(1986) gave examples of the growth of non-radial motions 
in N-body simulations. Arguments based on a numerical 
least-action method lead Peebles (1990) to the conclusion 
that irregularities in the mass distribution, together with 
external tides, induce non-radial motions that slow down 
the collapse. Lokas ct al. (1996) used N-body simulations 
and a weakly non-linear perturbative approach to study 
previrialization. They concluded that when the slope of 
the initial power spectrum is n > —1, non- linear tidal 
interactions slow down the growth of density fluctuations 
and the magnitude of the effect increases when n is 
increased. 

Opposite conclusions were obtained by Hoffman 
(1986a,1989), Evrard & Crone (1992), Bartelmann 
et al. (1993), Bcrtschinger & Jain (1994), Monaco 
(1995) and Kerscher et al. (2001). In particular Hoffman 
(1986a, 1989), using the quasi-linear (QL) approximation 
(Zcl'dovich 1970; Zcl'dovich & Novikov 1983) showed 
that the shear affects the dynamics of collapsing objects 
and it leads to infall velocities that are larger than in the 
case of non-shearing ones. Bcrtschinger & Jain (1994) 
put this result in theorem form, according to which 
spherical perturbations are the slowest in collapsing. 
Bartelmann ct al. (1993) argued that the collapse docs 
not start from a comoving motion of the perturbation, but 
that the continuity equation requires an initial velocity 
perturbation directly related to the density perturbation. 
The effect is that collapse proceeds faster than in the 
case where the initial velocity perturbation is set equal 
to zero and the collapse timcscalc is shortened. Kerscher 
et al. (2001) investigated the inhomogeneous, anisotropic 
collapse of non-isolated dust objects. Their collapse 
model, differently from the pure spherical collapse model, 
which focuses only on the mass overdensity, considers 
also the fluctuations in the expansion rate as well as 
the effect of shear fields on the collapse. Their results 
indicate that on large mass-scales the shear fields, caused 
by internal and external mass concentrations, accelerate 
the collapse. The N-body simulations by Evrard & Crone 
(1992) did not reproduce previriahzation effect, but the 
reason is due to the fact that they assumed an n = — 1 
spectrum, differently from the n = one used by Peebles 
(1990) that reproduced the effect. If n < — 1 the peculiar 
gravitational acceleration, g oc diverges at 

large R and the gravitational acceleration moves the fluid 
more or less uniformly, generating bulk flows rather than 
shearing motions. Therefore, its collapse will be similar 
to that of an isolated spherical clump. If n > —1, the 
dominant sources of acceleration are local, small-scale 
inhomogeneities and tidal effects will tend to generate 



non-radial motions and resist gravitational collapse. In 

a more recent paper, Audit et al. (1997) have proposed 
some analytic prescriptions to compute the collapse time 
along the second and the third principal axes of an ellip- 
soid, by means of the 'fuzzy' threshold approach. They 
pointed out that the formation of real virialized clumps 
must correspond to the third axis collapse and that the 
collapse along this axis is slowed down by the effect of 
the shear rather than be accelerated by it, in contrast 
to its effect on the first axis collapse. They concluded 
that spherical collapse is the fastest, in disagreement 
with Bcrtschinger & Jain's theorem. This result is in 
agreement with Peebles (1990). The quoted controversy 
was addressed by Del Popolo, Ercan & Xia (2001) who 
examined the evolution of non-spherical inhomogeneities 
in a Einstein-de Sitter universe, by numerically solving 
the equations of motion for the principal axes and the 
density of a dust ellipsoid. They showed that for lower 
values of u {u = 2) the growth rate enhancement of the 
density contrast induced by the shear is counterbalanced 
by the effect of angular momentum acquisition. For u > 3 
the effect of angular momentum and shear reduces, and 
the evolution of perturbations tends to follow the behavior 
obtained in the SM. Del Popolo (2002) studied the role 
of shear fields on the evolution of density perturbations 
by using an analytical approximate solution for the 
equations of motion of homogeneous ellipsoids embedded 
in a homogeneous background. The equations of motion 
of a homogeneous ellipsoid (Ickc 1973; White & Silk 1979 
(hereafter WS)) were modified in order to take account of 
the tidal field, as done in Watanabe (1993) and then were 
integrated analytically, similar to what was done in WS. 
The density contrast at turn-around and the collapse 
velocity were found to be reduced with respect to that 
found by means of the SM. The reduction increases with 
increasing strength of the external tidal field and with 
increasing initial asymmetry of the ellipsoids. 

The second physical effect with changes cluster col- 
lapse is dynamical friction. Former treatments of the dy- 
namical friction effects on the structure of clusters of 
galaxies, considering only the component generated by 
the galactic population on the motion of galaxies them- 
selves are due to White (1976) and Kashfinsky (1984, 
1986, 1987). AC recalculated the effect of dynamical fric- 
tion taking into account the effect of substructure, showing 
that dynamical friction delays the collapse of low-v peaks 
inducing a bias of dynamical nature. Because of dynam- 
ical friction under-dense regions in clusters (the clusters 
outskirts) accrete less mass with respect to that accreted 
in absence of this dissipative effect and as a consequence 
over-dense regions are biased toward higher mass (AC). 
Dynamical friction and non-radial motions acts in a simi- 
lar fashion: they delay the shell collapse consequently in- 
ducing a dynamical bias. Whenever efficient, these mech- 
anisms will generate a physical selection of those peaks in 
the initial density field that eventually will give rise to the 
observed cosmic structures. As a consequence of dynami- 
cal friction and tidal torques, one expects changes in the 
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threshold of collapse, the mass function and the correla- 
tion function. 

In this paper, I shall study how structure formation 
is changed by the joint effect of dynamical friction, tidal 
torques and a non-zero cosmological constant. The plan of 
the paper is as follows. In section 2, I summarize the role 
of dynamical friction in structure formation, in section 3 
the role of tidal torques. Section 4 shows how the joint 
effect of dynamical friction and tidal torques affect cluster 
formation and evolution. Section 5 deals with the modifi- 
cation of the threshold of collapse and section 6 with the 
modifications induced on the mass function. Section 7 is 
devoted to conclusions. 

2. Dynamical friction and structure formation. 

In a hierarchical structure formation model, the large 
scale cosmic environment can be represented as a collision- 
less medium made of a hierarchy of density fiuctuations 
whose mass, M, is given by the mass function N{M,z), 
where z is the redshift. In these models matter is con- 
centrated in lumps, and the lumps into groups and so 
on. In such a material system, gravitational field can be 
decomposed into an average field, Fo(r), generated from 
the smoothed out distribution of mass, and a stochastic 
component, Fstochir), generated from the fiuctuations in 
number of the field particles. The stochastic component of 
the gravitational field is specified assigning a probability 
density, W{F), (Chandrasekhar & von Neumann 1942). 
In an infinite homogeneous unclustered system VF(F) is 
given by Holtsmark distribution (Chandrasekhar & von 
Neumann 1942) while in inhomogeneous and clustered sys- 
tems iy(F) is given by Kandrup (1980) and Antonuccio- 
Delogu & Barandela (1992) respectively. The stochastic 
force, F stock, in a self-gravitating system modifies the mo- 
tion of particles as it is done by a frictional force. In fact 
a particle moving faster than its neighbors produces a de- 
fiection of their orbits in such a way that average den- 
sity is greater in the direction opposite to that of trav- 
eling causing a slowing down in its motion. Following 
Chandrasekhar & von Neumann's (1942) method, the fric- 
tional force which is experienced by a body of mass M 
(galaxy), moving through a homogeneous and isotropic 
distribution of lighter particles of mass m (substructure), 
having a velocity distribution n{v) is given by: 



M 



dt 



-AttG^ M^n{v)— log Ap 



(1) 



where log A is the Coulomb logarithm, p the density of the 
field particles (substructure). 

A more general formula is that given by Kandrup(1980) 
in the hypothesis that there are no correlations among 
random force and their derivatives: 



-77V 



JWiF)F^T{F)d^F 



2 < i;2 > 



(2) 



where r] is the coefficient of dynamical friction, T{F) the 
average duration of a random force impulse, < > the 



characteristic speed of a field particle having a distance 
r ~ (^^)^^^ from a test particle (galaxy). This formula 
is more general than Eq. ([I]) because the frictional force 
can be calculated also for inhomogeneous systems when 
W{F) is given. If the field particles are distributed homo- 
geneously the dynamical friction force is given by: 



F ■ 



-rjv 



< V 



log I 1.12 



(3) 



(Kandrup 1980), where ma and Ua are respectively the 
average mass and density of the field particles. Using virial 
theorem we also have: 



GmaUa 



Mt, 



(4) 



sys 



where Mtot is the total mass of the system, Rsys its radius 
and N is the total number of field particles. The dynamical 
friction force can be written as follows: 



F ■ 



-rjv 



AAAlGmaUac 



11/2 



where N = ^Rtysn 



J,/2 

^ is the comoving 



N 

and Uac — Ua X a 
number density of peaks of substructure of field particles. 
This last equation supposes that the field particles gener- 
ating the stochastic field are virialized. This is justified by 
the previrialization hypothesis (Davis & Peebles 1977). 
To calculate the dynamical evolution of the galactic com- 
ponent of the cluster it is necessary to calculate the num- 
ber and average mass of the field particles generating the 
stochastic field. 

The protocluster, before the ultimate collapse at z ~ 
0.02, is made of substructure having masses ranging from 
and from galaxies. I suppose that the stochas- 
tic gravitational field is generated from that portion of 
substructure having a central height v larger than a crit- 
ical threshold Vc- This latter quantity can be calculated 
(following AC) using the condition that the peak radius, 
fpkiv > Vc)i is much less than the average peak separa- 
tion Uai^ > i'c)~^^^, where Ua is given by the formula of 
BBKS for the upcrossing points: 



cxp(i/2/2) 7 



(27r) 



4^/3 



-^exp(-57^;.2/18)] 



(6) 



572(1-5779)1 

where 7, i?* are parameters related to moments of the 
power spectrum (BBKS Eq. 4.6A). The condition rpfc(i/ > 
Vc) < 0.1no(i^ > v^~^l^ ensures that the peaks of sub- 
structure are point like. Using the radius for a peak: 

r 1 1^^^ 



.(l + z/(To)(7^ + (0.9H)3/2_ 
(AC), I obtain a value of Vr = 1.3 and then we have 



v^) = 50.7Mpcr^ (7 = 0.4, R^ = bOKpc) and nia is given 
by: 

1 



nail" > Vc) 



(8) 
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(in accordance with the result of AC) , where rupk is given 
in Peacock & Heavens (1990) and Npk is the average 
number density of peak (BBKS Eq. 4.4). Galaxies and 
Clusters of galaxies are correlated systems whose autocor- 
relation function, ^(r), can be expressed, in a power law 
form (Peebles 1980; Bahcal & Soneira 1983; Postman et 
al. 1986; Davis & Peebles 1983; Gonzalez et al. 2002). The 
description of dynamical friction in these systems need to 
use a distribution of the stochastic forces, W{F), taking 
account of correlations. In this last case the coefficient of 
dynamical friction, 77, may be calculated using the equa- 
tion: 



T] = J d^FW{F)F^T{F)/{2 < v"^ >) 

and using Antonuccio & Atrio (1992) distribution: 
1 



W{F) 



2n^F 



dkksin{kF)Af{k) 



(9) 



(10) 



where Af, which is a linear integral function of the cor- 
relation function ^(r), is given in the quoted paper (Eq. 
36). 

3. Tidal torques and structure formation. 

The explanation of galaxies spins gain through tidal 
torques was pioneered by Hoyle (1949). Peebles (1969) 
performed the first detailed calculation of the acquisition 
of angular momentum in the early stages of protogalac- 
tic evolution. More recent analytic computations (White 
1984, Hoffman 1986, Ryden 1988a; Eisenstein & Loeb 
1995; Catelan & Theuns 1996a, b) and numerical simu- 
lations (Barnes & Efstathiou 1987) have re-investigated 
the role of tidal torques in originating galaxies angular 
momentum. One way to study the variation of angular 
momentum with radius in a galaxy is that followed by 
Ryden (1988a). In this approach the protogalaxy is di- 
vided into a series of mass shells and the torque on each 
mass shell is computed separately. The density profile of 
each proto-structure is approximated by the superposi- 
tion of a spherical profile, (5(r), and a random CDM dis- 
tribution, £(r), which provides the quadrupole moment of 
the protogalaxy. As shown by Ryden (1988a) the net rms 
torque on a mass shell centered on the origin of internal 
radius r and thickness 5r is given by: 

47r 



(11) 



(|rp)V2 ^ V30 ( -Gj [{a,m{rf){q^r^{rf) 
-(«2™(r)<?;„(r))2]i/2 



where g/m, the multipole moments of the shell and air, 
the tidal moments, are given by: 



(fen(r)2) = -^Mi / k'dkP{k)j,{krf 



{a2m(rY) = 



i2n) 



dkP (fc) ji {kry 



(12) 



(13) 



(a2™(r)g*„(r)) = -^PbM.u / kdkP {k) j^{kr) j2{kr){U) 



where Msh is the mass of the shell, ii{r) and j2{'r) are 
the spherical Bessel function of first and second order 
while the power spectrum P{k) is given by, BBKS (equa- 
tion (G3)): 



T{k) 



[In (1 + 2.Mq)] 



•[1 



2.34g 

-3.89g+ (16.1g)2 + (5.46g)3 + (6.71)''] ^^(15) 



(where q = n^^^^j— ^- Here 9 = pcr/(l-686p^) represents 
the ratio of the energy density in relativistic particles to 
that in photons {9 = 1 corresponds to photons and three 
flavors of relativistic neutrinos). The power spectrum was 
normalized to reproduce the observed abundance of rich 
cluster of galaxies (e.g., Bahcal & Fan 1998). Filtering 
the spectrum on cluster scales, Rf = 3h^^AIpc, I have 
obtained the rms torque, T(r), on a mass shell using Eq. 
(|lip then I obtained the total specific angular momentum, 
/i(r, J/), acquired during expansion integrating the torque 
over time (Ryden 1988a Eq. 35): 



/i(r, v) 



Toto -^-5/2 



sh 



(1 



cost 



1 

3 

.f2{^) 



2/3 



-dd 



(16) 



So 



the functions fi{^), f2{'&) are given by Ryden (1988a) (Eq. 
31) while the mean over-density inside the shell, 6{r), is 
given by Ryden (1988a): 



6{r, v) 



daa^6{a) 



(17) 



where S{r) 



As showed by Ryden (1988a), the 



P(r)-Pb 
Pb 

rms specific angular momentum, h{r,v)^ increases with 
distance r while peaks of greater u acquire less angular mo- 
mentum via tidal torques. This is the angular momentum- 
density anticorrelation showed by Hoffman (1986). This 
effect arises because the angular momentum is propor- 
tional to the gain at turn around time, tm, which in turn 
is proportional to (5(r, i/)^^ cx u^^/'^ . 



4. IVIodification of collapse. 

Tidal torques and dynamical friction acts in a similar fash- 
ion. As reported in the introduction AC calculated the ef- 
fect of dynamical friction taking into account the effect of 
substructure, showing that dynamical friction delays the 
collapse of low-i^ peaks inducing a bias of dynamical na- 
ture. Similarly non-radial motions would slow the rate of 
growth of the density contrast by lowering the peculiar 
velocity and suppress collapse once the system detaches 
from general expansion. In fact, in the central regions of 
a density peak [r < 0.5Rf) the velocity dispersion at- 
tain nearly the same value while at larger radii (r > Rf) 
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remembering that 




0.2 0.4 



Fig. 1. The time evolution of the expansion parameter. 
The solid hne is a(t) for the SM; the dashed line is a(t) 
taking account only dynamical friction; the dotted line is 
a(t) taking account of the cumulative effect of non-radial 
motions and dynamical friction in the case of a, v = 2 
peak. 

the radial component is lower than the tangential compo- 
nent. This means that motions in the outer regions are 
predominantly non-radial and in these regions the fate of 
the infalling material could be influenced by the amount 
of tangential velocity relative to the radial one. This can 
be shown writing the equation of motion of a spherically 
symmetric mass distribution with density n(r) (Peebles 
1993): 



d_ 

dt 



0(18) 



where (vr) and (v^) are, respectively, the mean radial and 
tangential streaming velocity. Eq. ^TE\i shows that high 
tangential velocity dispersion ((w^) > 2{vf)) may alter 
the infall pattern. The expected delay in the collapse of a 
perturbation, due to non-radial motions, dynamical fric- 
tion and also taking account of a non-zero cosmological 
constant, may be calculated solving the equation for the 
radial acceleration (Kashlinsky 1986, 1987; AC; Peebles 
1993): 



dVr 
IF 



L^jr, v) 



. . dr A 



(19) 



where L{r, v) is the angular momentum, g{r) the accelera- 
tion, and A the cosmological constant. Writing the proper 
radius of a shell in terms of the expansion parameter, 
a{ri,t): 



M 



(21) 



and that pf, = where Hq is the Hubble constant and 
assuming that no shell crossing occurs so that the total 
mass inside each shell remains constant, that is: 



the Eq. may be written as: 



dt2 



2a2 



(22) 



da A 
^Tt + 3 « 



Eq. does not take into account back-reaction terms, 
which have been taken into account by Kerscher et al. 
(2001). On large mass-scales this generic model leads to 
changes in the pure top-hat model. One consequence of 
the Kerscher et al. (2001) model is the prediction of more 
collapsed objects at high-mass scales, which, as we shall 
see in the next sections, leads to conclusions similar to the 
model of the present paper. 

The equation (j23p may be solved using the initial con- 
ditions: (^) = 0, a = ttraax — 1/(5 and using the function 
/i(r, v) = L{r, v)/Msh found in Section 2 to obtain a(t) 
and the time of collapse, Tc{r, u). As shown by Gunn & 
Gott (1972), this last quantity in the case of a pure SM 
(namely when tidal torques, dynamical friction and cos- 
mological constant are not taken into account) is given 
by: 



(24) 



r{ri,t) = ria{ri,t) 



(20) 



In Fig. 1, I show the effects of non-radial motions and dy- 
namical friction separately, in the case of a = 2 peak. 
As displayed non-radial motions have a larger effect on the 
collapse delay with respect to dynamical friction. In Fig. 
2, I compare the results for the time of collapse, Tc, for 
= 2, 4 with the time of collapse of the classical SM (Eq. 
21]) . As shown the presence of non-radial motions produces 
an increase in the time of collapse of a spherical shell. The 
collapse delay is larger for low value of v and becomes 
negligible for > 3. This result is in agreement with the 
angular momentum-density anticorrelation effect: density 
peaks having low value of v acquire a larger angular mo- 
mentum than high v peaks and consequently the collapse 
is more delayed with respect to high v peaks. 

Given Tc{r, v), I also calculated the total mass gravita- 
tionally bound to the final non-linear configuration. There 
are at least two criteria to establish the bound region to 
a perturbation 6{r): a statistical one (Ryden 1988b), and 
a dynamical one (Hoffman & Shaham 1985), summarized 
in the following. 

In biased galaxy formation theory structures form 
around the local maxima of the density field. Every den- 
sity peak binds a mass M that can be calculated when we 
know the binding radius of the density peak. The radius of 
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Fig. 2. Fig. 2a. The time of collapse of a shell of matter 
in units of the age of the universe to ior v — 2 (dotted 
line) compared with SM (solid line). Fig. 2b. The time 
of collapse of a shell of matter in units of the age of the 
universe to for = 4 (dotted line) compared with SM 
(solid line). 



the bound region for a chosen density profile 5{r) may be 
obtained in several ways. A first criterion is statistic. The 
binding radius of the region, r;,, is given by the solution 
of the equation: 



< 5{r) >=< {6- <d>f > 



1/2 



(Ryden 1988). At radius r << ri, the motion of particles 
is predominant toward the peak while when r >> ri, the 
particle is not bound to the peak. Another criterion that 
can be used is dynamical. It supposes that the binding 
radius is given by the condition that a shell collapse in a 
time, Tc, smaller than the age of the universe to- 



Tc{r) < to 



(26) 



(Hoffmann & Shaham 1985). This last criterion, differ- 
ently from the previous one, contains some prescriptions 
particularly connected with the physics of the collapse pro- 
cess of a shell. For this reason I used it to calculate the 
binding radius. The time of collapse, Tc{r), at radius r 
can be obtained solving numerically Eq. for differ- 
ent values of Si, the initial overdensity, from a given den- 
sity profile 6{r). I use the average density profile given by 
BBKS: 



S{r) ^ A 



[ <(r) 

\C(0)l/2 ^^(0)1/2(1-^2) 



6'(z/7,7) 



where A is a constant given by the normalization of the 



perturbation spectrum, P{k), 



a(M) 



where Sm — 



(25) 



1.686 is the critical threshold for a SM, a{M) is the r.m.s. 
density fluctuation on the mass scale M, ^(r) is the cor- 
relation function of two points, 7 and two constants 
obtainable from the spectrum (see BBKS) and finally 
0{'yi>, 7) is a function given in the quoted paper (eq. 6.14). 
Given the average density profile the average density in- 
side the radius r in a spherical perturbation is given by 
Eq. ini). 

Finally, I calculated the binding radius, rb{i'), for a 
SM, calculating rco(r-) by means of Eq. p4|) and the den- 
sity profile given in Eq. (|27|) and then applying the con- 
dition Tco(r) < to- I repeated the calculation for 1.7 < 
I' < 4. Then I repeated the calculation using T'c(r), the 
collapse time that takes into account non-radial motions 
and dynamical friction. I found a relation between v and 
the mass of the cluster using the equation: M — ^r^pi,. 
The result is the plot in Fig. 3 for the binding radius r;, 
versus v. 

In fig. 3, I compare the peak mass obtained from SM, 
using Hoffman & Shaham's (1985) criterion, with that ob- 
tained from the model taking into account non-radial mo- 
tions, dynamical friction and A 0. As shown for high 
values oiv {v > 3) the two models give the same result for 
the mass while for v < i the effect of non-radial motions 
produces less bound mass with respect to SM. Decreasing 
the effect of non-radial motions produces a decrease in the 
bound mass. 

The situation represented in the previous three figures 
may be summarized as follows: dynamical friction and 
non-radial motions delays the collapse of perturbations. 
Both effects act in the direction of delaying structure col- 
lapse, so that their effects add. The effects have a similar 
magnitude, but non-radial motions induce a slight larger 
delay in collapse. As a consequence of this delay of collapse 
the matter bound to structures is less than what expected 
in the case of SM. 
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5. The threshold of collapse Sc- 

In this section, I am going to show how dynamical friction 
and tidal fields influence the critical overdensity threshold 
for the collapse, Sc, which is not constant as in a SM but it 
depends on mass. An analytic determination of (5c (i^) can 
be obtained following a technique similar to that used by 
Bartlett & Silk (1993). 



2.6 - 



2.4 - 




1 1.5 2 2.5 



Fig. 3. Variation of the binding radius rf, with i^. The solid 
line is the binding radius in the SM, while the dashed line 
is the same as in presence of non-radial motions dynamical 
friction and A 7^ 0, for = 2. 

Using Eq. it is possible to obtain the value of 

the expansion parameter of the turn around epoch, amax , 
which is characterized by the condition ^ = 0. Using the 
relation between v and Si, in linear theory (Peebles 1980), 
I find: 

BiM) = 5, = Sco 

where Sco = 1-686 is the critical threshold for SM, n is the 
initial radius, rta is the turn-around radius, Xo ~ e-oTco and 
ix{5) is given in Colafrancesco, Antonuccio & Del Popolo 
(1995) (Eq. 29). The quantity L appearing in Eq. ^ 
is the total angular momentum acquired by the proto- 
structure during evolution. In order to calculate L, it is 
possible to use the same model described in Sec. 3. (more 
hints on the model and some of the model limits can be 
found in Del Popolo, Ercan & Xia 2001). The result of the 
calculation is shown in Fig. 5, where I plot Sc{v) obtained 
by means of the model of the present paper together with 



that obtained by ST using an ellipsoidal collapse model. 
The dashed line represents Sc{v) obtained with the present 
model, while the solid line that of ST. Both models show 
that the threshold for collapse decreases with mass and 
when > 3 the threshold assume the typical value of the 
SM. In other words, this means that, in order to form 
structure, more massive peaks must cross a lower thresh- 
old, 5c{v), with respect to under-dense ones. At the same 
time, since the probability to find high peaks is larger in 
more dense regions, this means that, statistically, in order 
to form structure, peaks in more dense regions may have 
a lower value of the threshold, Sc{v), with respect to those 
of under-dense regions. This is due to the fact that less 
massive objects are more influenced by external tides, and 
consequently they must be more overdense to collapse by 
a given time. In fact, the angular momentum acquired by 
a shell centered on a peak in the CDM density distribution 
is anti-correlated with density: high-density peaks acquire 
less angular momentum than low-density peaks (Hoffman 
1986; Ryden 1988). A larger amount of angular momen- 
tum acquired by low-density peaks (with respect to the 
high-density ones) implies that these peaks can more eas- 
ily resist gravitational collapse and consequently it is more 
difficult for them to form structure. This is in agreement 
with Audit et al. (1997), Peebles (1990), which pointed 
out that the gravitational collapse is slowed down by the 
effect of the shear rather than fastened by it (as sustained 
by other authors). Therefore, on small scales, where the 
shear is statistically greater, structures need, on average, 
a higher density contrast to collapse. This results in a 
tendency for less dense regions to accrete less mass, with 
respect to a classical SM, inducing a biasing of over-dense 
regions towards higher mass. 

6. Mass function 

In cosmology, it is possible to model the statistics of halo 
formation and mergers, in a theory sometimes referred 
to as extended-Press-Schechter, or excursion set (see PS, 
Bond et al. 1991, Lacey & Cole 1993) by means of a ran- 
dom walk. As well known, a random walk is a stochastic 
process consisting of a sequence of uncorrelated discrete 
steps. The idea can be summarized as follows. Suppose 
that density perturbations are represented by a density 
field and imagine smoothing the density field around this 
point with a filter, say a top-hat filter, of a successively 
smaller radius. In a hierarchical universe such as our own, 
one expects that a plot of the smoothed density versus 
the size of the filter traces out a random walk. The overall 
trend is for the smoothed density to rise with a smaller 
filter radius, but with significant fluctuations. 

The excursion set theory postulates a flat threshold 
or barrier (usually motivated by spherical collapse): once 
the random walk crosses the barrier, a halo is declared 
to have formed with a size given by the crossing radius. 
Often the filter radius R is denoted by alternative equiv- 
alent quantities, e.g.: the mass variance S. The latter is 
particularly useful for Gaussian random initial conditions 
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Fig. 4. The mass accreted by a collapsed perturbation, in 
units of 10^^ Mq, taking into account non-radial motions, 
dynamical friction and a non-zero cosmological constant 
(dotted line) compared to SM mass (solid line). 



- note that S increases as the radius or mass is decreased. 
A natural question to ask is: given Gaussian random ini- 
tial conditions, what is the probability that a random walk 
first crosses the barrier between S and S + dS, denoted by 
f{S)dSl The quantity f{S) is known as the first-crossing 
distribution. In excursion set theory, the first-crossing dis- 
tribution is directly related to the halo mass function, a 
key quantity of interest. 

Some years ago, it was found that it useful to consider 
random walks with a non-flat barrier, which is usually re- 
ferred to as moving barrier, in the sense that the barrier 
'moves' with S. For instance, Sheth et al. (2001, SMTOl 
hereafter) found that an ellipsoidal collapse model sug- 
gests a moving barrier for defining halos, which produces 
a first-crossing distribution, or halo mass function, that 
better matches N-body simulations. SMTOl estimated the 
mass function assuming an anisotropic ellipsoidal collapse, 
where the time evolution of the half-axes is determined by 
the Zeldovich approximation. Kerscher et al. (2001) gen- 
eralized this ansatz to domains not restricted to be ellip- 
soidal throughout the evolution. Several other approaches 
are also based on the Zeldovich approximation, sometimes 
the third-order Lagrangian perturbation approximation is 
employed. In these local approaches the collapse of a do- 
main is associated with a diverging local matter density 
(see e.g. Bartelmann et al. 1993, Monaco 1997a, b, Lee & 
Shandarin 1998.) All these models, similarly to the re- 
sults obtained in the present paper, for the mass func- 
tion, agrees on predicting a larger abundance of collapsed 
objects at high-mass scales, similarly to N-body simula- 



Fig. 5. The critical threshold, (5c(!^) versus v. The dashed 
line is obtained with the model of the present paper while 
the solid line is that of ST. 

tions like those of Jenkins et al. 2001 (JOl); R03; Yahagi, 
Nagashima & Yoshii (2004) (YNY); Warren et al. 2005 
(W05)). STl connected the form of the barrier with the 
form of the mass function. As shown by STl, for a given 
barrier shape, B{S), the first crossing distribution is well 
approximated by: 



/(5).5HT(S)|»p,-a|)!)|| 



(29) 



where T{S) is the sum of the first few terms in the Taylor 
expansion of B{S): 



T{S) 



n=0 



(30) 



In the case of the barrier shape given in Eq. ([28]) the Eqs. 

give, after truncating the expansion at n = 5 
(see ST), the multiplicity function, i'f{i'): 

/3ig{ai) ^25(^2) , [339(^3) 



yf{y) = ( 1 
exp{- 



[avf 



I3l , P2 , /33 




where 



g{ai) =1 l~ai+ 



ai{ai 



1 



(aj/)"^ [au)"^ {avT' 
) a,(a, - 1) • • • (a, - 4) 



2! 5! 

where i = 1 or 2, ai = 0.585, j3i = 0.46, a2 — 0.5 and 
P2 = 0.35, a3 = 0.4 and /Sg = 0.02, and a = 0.707. 

The "multiplicity function" is correlated with the 
usual, more straightforwardly used, "mass function" as 



(32) 
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follows. Following ST notation, if f(M,S)dM denotes the 
fraction of mass that is contained in collapsed haloes that 
have mass in the range M-M + dM , at redshift z, and S{z) 
the redshift dependent overdensity, the associated "uncon- 
ditional" mass function is: 



n{M, S)dM = j^f{M, 6)dM 



(33) 



In the excursion set approach, the universal or "uncon- 
ditional" mass function, n{M, z), representing the average 
comoving number density of haloes of mass M is given by: 



n(M, z) 



Pb{z) dXogv 
M2 dlogM 



(34) 



(Bond et al. 1991). The function vf{v) is the "multiplic- 
ity function" which is obtained by computing the distri- 
bution of first crossings, f{v)dv, of a barrier B{v), by in- 
dependent, uncorrelated Brownian motion random walks. 
Multiplicity function and mass function are related by Eq. 
(|34|) . It is to be noted that in literature sometime the terms 
mass function and multiplicity function are used as syn- 
onymous (e.g. ST, Lin et al. 2002). 

With a similar calculation, STl found, that 



except at very high redshifts, where it significantly over- 
predicts the halo abundance. At all redshifts up to z = 10, 
the difference is < 10% for each of our well sampled mass 
bins. However, the ST function begins to overpredict the 
number of haloes increasingly with redshift for z>10, up 
to ~ 50% by z=15. The simulation mass functions ap- 
pear to be generally steeper than the ST function, espe- 
cially at high redshifts. The ST function fits the simulated 
mass function to better than 10% over the range of -1.7 
< luCT"! < 0.5 while it appears to significantly overpre- 
dict haloes for h\a~^ > 0.5. JOl also found an overpre- 
diction by the ST function for Intr"^ > 0.75, which with 
their larger simulation volumes, corresponded primarily 
to objects of z<2 and of much higher mass. Additionally, 
JOl found the mass function to be invariant with redshift 
within their own results. 

R03 also considered the possibility of an empirical ad- 
justment to the ST function. They inserted a crude mul- 
tiplicative factor to the ST function as follows, with 5c = 
1.686 and FOE U =0.2: 



/(a) = /(a;ST) 



exp[-0.7/{a[cosh{2a)]^)] 



(39) 




exp{- 



0.5 



[av) 



0.6 



with ^2 — 1- This last result is in good agreement with 
the fit of the simulated first crossing distribution (ST): 



vf{v)dy - A3 1 + 




[auf 



exp(— aj//2) 



(36) 



where p = 0.3, and a = 0.707. 

The normalization factor ^3 has to satisfy the con- 
straint: 



valid over the range of -1.7 < Incr^^ < 0.9. The resulting 
/2fji^i36l)ion is virtually identical to the ST function for all 
—00 < Incr < 0.4. At higher values of Incr^^, this function 
declines relative to the ST function, reflecting an under- 
abundance of haloes that becomes greater with increasing 
Incr-i. For -1.7 < Ina^^ < 0.5, Eq. ^ matches R03 
data to better than 10% for well-sampled bins, while for 
0.5 < Incr"^ < 0.9, where poisson errors are larger, data 
is matched to roughly 20%. Note that in their notation, 
similarly to JOl: 



M dn{M,z) 



f{y)dv = 1 



(37) 



and as a consequence it is not an independent parameter, 
but is expressed in the form: 



Pb{z) diner 1 
Using the relation v ~ (%^)^ and Eq. 
f{a,z) = 2yf{y) 



(40) 

one finds that: 
(41) 



A. 



1 + 2-^t:-^''^T{1/2-p) 



0.3222 



(38) 



Recently in order to investigate the functional form of 
the universal multiplicity function, several authors (JOl; 
R03; YNY; W05) performed runs of N-body simulations 
with high mass resolution and compared them with differ- 
ent multiplicity function. The result of the previous stud- 
ies is that even the ST multiplicity function does not give 
correct prediction in some ranges of mass. For example, 
R03 used a high resolution ACDM numerical simulation 
to calculate the mass function of dark matter haloes, and 
its evolution, down to the scale of dwarf galaxies, back to 
a redshift of fifteen, in a 50 /i~^Mpc volume containing 
80 million particles. Their low redshift results allow us to 
probe low a density fiuctuations significantly beyond the 
range of previous cosmological simulations. They showed 
that the ST function provides a good fit to their data. 



The evolution of the mass function can be calculated 
evaluating it at different redshifts. The mass function, Eq. 

depends on redshift because of the dependence of 
Pb(z), V or a from z. 

For example cr(Af, z) = (t(M, z = 0)&(z), where h{z) 
evolves as (l+z)^^ in an f^o = 1 universe, and more slowly 
in a ACDM universe. In the following I shall plot the mass 
function of the present paper and I compare it with R03 
results. 

In Fig. 6, 1 plot the mass function for all of our outputs 
in the f{a) — ln(CT~^) plane. Large values of Incr"^ corre- 
spond to rare haloes of high redshift and/or high mass, 
while small values of lncr~^ describe haloes of low mass 
and redshift combinations. The solid line is the ST mass 
function while the dashed line the one obtained in the 
present paper and the dotted line represents Eq. (|39p. The 
ST function fits the simulated mass function to better than 
10% over the range of -1.7 < Ina^^ < 0.5 while it appears 
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Fig. 6. Mass function plotted in redshift independent form 
for all of R03 outputs: redshifts used are 0, 1., 2., 3., 4., 
5., 6.2, 7.8, 10., 12.1, 14.5. The solid line is ST prediction 
while the dashed and dotted line represent the result of 
the present paper and Eq. p9|) . respectively. 

to significantly overpredict haloes for In cr"^ > 0.5. The 
magnitude of the ST overprediction at high values of Ina^^ 
is consistent with being a function purely of Incr^^ rather 
than redshift, a natural consequence of the fact that the 
mass function is self similar in time (e.g Efstathiou et al. 
1988; Lacey & Cole 1994; JOl). JOl also found an overpre- 
diction by the ST function for Ina^^ > 0.75, which with 
their larger simulation volumes, corresponded primarily 
to objects of z<2 and of much higher mass. Additionally, 
JOl found the mass function to be invariant with redshift 
within their own results. The empirical adjustment to the 
ST mass function (Eq. ([55]) ). dotted line, describes much 
better numerical simulations data: for -1.7 < In cr^^ < 0.5, 
Eq. (I39p matches R03 data to better than 10% for well- 
sampled bins, while for 0.5 < Incr"^ < 0.9, where poisson 
errors are larger, data is matched to roughly 20%. The ST 
and the mass function of the present paper differs more 
in the high mass region, where the mass function of the 
present paper is steeper than ST and in better agreement 
with numerical simulations data than ST mass function. 
The better agreement between R03 simulations and the 
result of the present paper is connected to the improved 
barrier taking explicitly account the effects of several phys- 
ical effects: namely dynamical friction, tidal torques and 
the cosmological constant. 

7. Conclusions 

In this paper, I have studied the role of non-radial motions, 
dynamical friction and a non-zero cosmological constant 



on the collapse of density peaks solving numerically the 
equations of motion of a shell of baryonic matter falling 
into the central regions of a cluster of galaxies. I have 
shown that non-radial motions and dynamical friction pro- 
duce a delay in the collapse of density peaks having low 
value of v while the collapse of density peaks having > 3 
is not influenced. Both effects produce a delay in collapse, 
but the effect of non-radial motions is slightly larger than 
that of dynamical friction. A first consequence of this ef- 
fect is a reduction of the mass bound to collapsed pertur- 
bations and a raising of the critical threshold, 5c, which 
now is larger than that of the top-hat spherical model 
and depends on v. This means that shells of matter of 
low density have to be subjected to a larger gravitational 
potential, with respect to the homogeneous SM, in or- 
der to collapse. The delay in the proto-structures collapse 
gives rise to a dynamical bias similar to that described in 
CAD. Finally, I used the Extended Press-Schechter (EPS) 
and STl approach to calculate the mass function starting 
from the threshold function. Then, I compared the numer- 
ical mass function given in R03 with the theoretical mass 
function obtained by means of the excursion set model and 
an improved version of the barrier shape obtained in the 
present paper, which implicitly takes account of dynamical 
friction, tidal interactions between clusters and a non-zero 
cosmological constant. I showed that the barrier obtained 
in the present paper gives rise to a better description of 
the mass function evolution with respect to other models 
(ST, STl) and that the agreement is based on sound the- 
oretical models and not on fitting to simulations like in 
R03, YNY, or Warren et al. (2005). 

The main results of this last part of the paper can be 
summarized as follows: 

1) the non-constant barrier of the present paper combined 
with the STl model gives a mass function evolution in 
better agreement with the N-body simulations of R03 
than other previous models (ST, STl). 

2) The mass function of the present paper gives a good 
fit to simulations results as the fit function proposed by 
R03, but differently from that it was obtained from a 
sound theoretical background. 

3) The excursion set model with a moving barrier is very 
versatile since it is very easy to introduce easily several 
physical effects in the calculation of the mass function 
and its evolution, just modifying the barrier. 

4) The behavior of the mass function and its evolution at 
small masses is similar to that of ST, STl, but at higher 
values of mass or redshift it is steeper than ST, STl in 
agreement with N-body simulations of R03. 

The above considerations show that the excursion set 
approach that incorporates a non-spherical collapse which 
takes account of angular momentum acquisition, dynam- 
ical friction and a non-zero cosmological constant gives 
accurate predictions for a number of statistical quantities 
associated with the formation and clustering of dark mat- 
ter haloes. The improvement is probably connected also to 
the fact that incorporating the non-spherical collapse with 
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increasing barrier in the excursion set approach results in 
a model in which fragmentation and mergers may occur, 
effects important in structure formation. Moreover, the ef- 
fect of a non-zero cosmological constant adds to that of 
angular momentum and dynamical friction slightly chang- 
ing the evolution of the multiplicity function with respect 
to open models with the same value of matter density pa- 
rameter. 
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